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We explore the response of an unconventional superconductor to spatially inhomogeneous 
antiferromagnetism (SIAFM). Symmetry allows the superconducting order parameter in the E- 
representation models for UPt3 to couple directly to the AFM order parameter. The Ginzburg- 
Landau equations for coupled superconductivity and SIAFM are solved numerically for two possible 
SIAFM configurations: (I) abutting antiferromagnetic domains of uniform size, and (II) quenched 
random disorder of 'nanodomains' in a uniform AFM background. We discuss the contributions 
to the free energy, specific heat, and order parameter for these models. Neither model provides a 
satisfactory account of experiment, but results from the two models differ significantly. Our results 
demonstrate that the response of an E2 U superconductor to SIAFM is strongly dependent on the 
spatial dependence of AFM order; no conclusion can be drawn regarding the compatibility of E2 U 
superconductivity with UPt3 that is independent of assumptions on the spatial dependence of AFM. 
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I. INTRODUCTION 



The nature of the spatially inhomogeneous small- 
moment antiferromagnatism (AFM) observed in neutron 
scattering experimentsElTJ below T/v ~ 6K and its in- 
teraction with superconductivity (T c ~ 0.5K) remain 
central issues in determining the symmetry of the su- 
perconducting order- .parameter of UPt3. The unusual 
H-T phase diagramot] at ambient pressure apparently 
shows three superconducting phases in the mixed state 
and two Meissner phases. Experimental studies using 
hydrostaticO and uniaxials pressure reveal the existence 
of a critical pressure above which the zero-field transi- 
tion splitting disappears. This complex phase diagram 
strongly suggests that superconductivity in this heavy- 




has provided mo- 
roposed theories 



electron material is unconventiona 
tivation for much theoretical work. 

range from odd-in-frequency pairingiil to multicompo 
nent order parameters. The latter may belong to a single 
multidimensional representation of the symmetry group 
(see Rcfs. [?],|]for reviews), or they may belong to different 
representations of the crystal point group that are either 
accidentalU^degeneratejij reflect a 'higher symmetry' of 
the crystafcj or of spin spaceo. 

As briefly summarized below, experiments suggest an 
intriguing interplay between superconductivity and anti- 
ferromagnetism; however, an interpretation without sig- 
nificant ambiguity has not yet emerged. Here, wc focus 
on one theoretical proposal: an odd-parity superconduct- 
ing state with a two-dimensional order parameter that 
transforms like a single representation, the E 2 represen- 



tation of the hexagonal symmetry group. This order pa- 
rameter may also couple to the AFM order parameter. 
While the spatially homogeneous superconducting states 
of this model and their response to an applied magnetic 
field have been studied, there is comparatively little work 
that explores the effect of spatially inhomogeneous anti- 
ferromagnetism (SIAFM) on superconductivity. Moti- 
vated by recent work of Garget, which concludes that 
spatially varying AFM would rule out the two dimen- 
sional representation models for UPt3, we examine the 
sensitivity of an E 2u superconductor to SIAFM through 
numerical calculations of its response to two qualitatively 
different kinds of spatial configurations of the SIAFM. 
Since the free energy functionals for an E2« and an Ei g 
superconductor are formally identical, our results are also 
relevant to Ei g superconductors. 

Antiferromagnetism is the prime suspect for inducing 
the zero-field double phase transition observed in spe- 
cific heat experiments on high-quality crystals. A cou- 
pling between AFM and superconducting order parame- 
ters is suggested by a downward kink in the magnitude of 
the modulus of the AEM order parameter at the super- 
conducting transitionn In a comparison of specific heat 
measurements with neutron scattering experiments un- 
der pressure, Hayden et al.u observed that the disappear- 
ance of the double transition is correlated with the dis- 
appearance of signatures of antiferromagnetism. Recent 
work of Keizcr et alx-i further supports this correlation. 
Upon substituting Pd for Pt on a small number of sites, 
they find that the magnetic moment and the splitting 
of the double transition increase simultaneously with in- 
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creasing Pd doping. 

The nature of the antifcrromagnctism is itself unusual. 
A signature of an AFM phase transition in thermody- 
namic, NMR (nuclear magnetic resonance), and zero-field 
/iSR (muon spin, relaxation) experiments has so far not 
been observed.Ej This has been taken as evidence for the 
absence of long-range order and the existence of mag- 
netic fluctuations on a characteristic scale smaller than 
neutron scattering frequencies, but greater than those 
of NMR. This temporal fluctuation has largely been ig- 
nored and the antiferromagnetism has been taken to be 
static when considering the interaction of AFM with 
superconductivity. While conventional thermodynamic 
signatures of the Neel temperature, TV, have not been 
observed, more recent transverse high-field //SR exper- 
iments have detected anomalies at TV as identified by 
neutron scattering.Ej 

The appearance of the double transition was unex=. 
pected. As shown in Fig. |l|, specific heat experimentstZl 
prior to 1989 typically showed a single anomalously broad 
peak at the transition to superconductivity. An obvious 
explanation is that AFM domains increase in size during 
the annealing process sharpening the distribution around 
two intrinsic superconducting transitions, but X-raya and 
neutron scatteringta experiments fail to show any obvi- 
ous correlation between domain size and annealing. 

A commonly held physical picture is based on an inter- 
pretation of the neutron scattering data. The AFM or- 
der with orthorhombic symmetry appears with ~ 0.02/is 
ordered moments constrained to lie in the basal plane. 
AFM order occurs in domains of uniform size ~ 30 nm, 
that are randomly distributed over three q vectors that 
are oriented at 0°, 120°, and 240° with respect to the a* 
axis. The moments are essentially rigidly locked to the 
lattice for fields in excess of the zero-tempatature upper 
critical field of the superconducting state.t3E3 

This picture is not firmly established. Existing neutron 
scattering data is unable to rule out the possibility that 
instead of domains of a singlc-q structure, AFM order ap- 
pears in a triplc-q stpjeture that preserves the symmetry 
of the crystal lattice. EHl Moreover;, a recent careful analysis 
of the neutron scattering dataO finds that no conclusion 
can be drawn from existing data on whether the stag- 
gered moment remains fixed to the lattice or whether it 
rotates with an applied magnetic field. Existing data also 
cannot distinguish randomly oriented abutting domains 
with small magnetic moments from small domains with 
magnetic moments of ~ 1[j,b interspersed in an otherwise 
nonmagnetic system. 

Whether AFM occurs in a few large domains with stag- 
gered moments that are free to rotate in an applied field, 
or whether it is strongly spatially varying with staggered 
moments that are rigidly fixed to the lattice, is impor- 
tant in determining the symmetry of the superconduct- 
ing states, particularly for states near the upper tran- 
sition temperature. Several authors take the view that 
this quasistatic AFM acts as a symmetry breaking field 
(SBF) and lifts the degeneracy among components of a 



multicomponcnt order parameter resulting in two super- 
conducting phase transitions acpacated by ~ 50 mK that 
are observed in zero field BE3~EJQ'E3i23 A two-component 
odd-parity order parameter that transforms like the two- 
dimensional Eiu representation of the hexagonal symme=j 
try group Deh is one of the more promising proposals .Q'EU 
At low temperature or in the absence of the SBF, weak 
coupling BCS theory shows that the homogeneous equi- 
librium state breaks time-reversal symmetry. Further 
calculations using weak-coupling BCS theory show that 
thermal conductivity,^ -transverse sound attenuation,Ej 
and upper critical fieldsE-3 of this state are in good agree- 
ment with experiments for temperatures in the low tem- 
perature phase£3 Coupling to an SBF has been included 
within a Ginzburg-Landau (GL) theory developed for 
a single-domain superconducting state. Signatures of a 
double phase transition-ate apparent in the specific heat 
and lower critical field and in the cores of vortices SB 
In contrast to two-dimensional even-parity Ei s and E2 9 , 
and to odd-parity Ei„ order parameters, the E2„ model 
can allow for a tetracritical point for arbitrary field ori- 
entations in the H-T phase diagram. An enhancement 
of this model includes the competition between magnetic 
anisotropyr-apd Zccman energies of the magnetic order 
parameter J23 and reproduces the angular dependence of 
the upper critical field observed in experiment E3 

Within the E2 U model, comparatively little has been 
done to explore the effect of coupling superconductiv- 
ity and spatially inhomogeneous AFM. Motivated by the 
quenched domain interpretation of the neutron scattering 
data, early wock. on E-representation superconductivity 
by Joynt et al!El and by MineevEa focused on abutting 
AFM domains that are uniform in size with dimensions 
of a superconducting coherence length and considers the 
possibility of a superconducting glass phase. Based on 
a variational calculation for domains of uniform size and. 
on calculations for a one-dimensional 'toy model', GargE3 
argued that the pure E-representation models are 'incom- 
patible' with UPt3 for the small domains suggested by 
the neutron scattering experiments. 

Taking the AFM order to be static, we explore the 
sensitivity of E2„ superconductivity to spatially varying 
AFM order in two models for the (disordered) domain 
structure of the AFM state: (I) abutting AFM domains 
of uniform size with orientations that are equally dis- 
tributed over the allowed q-vectors, and (II) small do- 
mains interspersed in a homogeneous AFM background. 
We consider these models to represent limiting cases for 
the configuration of the SIAFM. The first model cor- 
responds to the standard interpretation of the neutron 
scattering data. The second model begins with uniform 
AFM order and adds 'nanodomains' with random posi- 
tion and orientation of the staggered magnetization; this 
provides the broadening mechanism for the linewidth of 
neutron scattering data. We present numerical solutions 
of the GL equations in two spatial dimensions and fo- 
cus on contributions to the free energy, the specific heat, 
and the nature of the superconducting state. Our results 
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for Model I agree with those for Garg's simple periodic 
model. While neither model provides an adequate ac- 
count of experiment, our results for Models I and II dif- 
fer, and taken together they do not lead to Garg's con- 
clusion that the E% u model is incompatible with UPt3. 
Rather, they suggest that a reasonably accurate model 
of the spatial dependence of the AFM is required to make 
meaningful comparisons with experiment. In a forthcom- 
ing work, we discuss other domain configurations along 
with the effects of dimensionality, and the possible role 
of superconducting glass phases.c3 Central results of this 
paper are contained in the comparison of specific heat 
calculations with experiment for both models. 

In Model I, the signatures of the superconducting tran- 
sitions rapidly smear and broaden with decreasing do- 
main size. The SBF introduces a convenient length scale 
£ e (defined below) which, for UPt3, is some three times 
larger than the zero-temperature superconducting co- 
herence length. For domain sizes of 10£ e — 20£ e , the 
calculated specific heat compares well with data from 
experiment Ji3 This homogeneous domain size model does 
not agree with neutron scattering experiments which, 
when viewed through the lens of Model I, would suggest 
a much smaller domain size, ~ l£ e — 3£ e . As the domain 
size is decreased below ~ 2 £ £ only a single supercon- 
ducting transition appears in the specific heat. While 
the suppression of the double transition with decreas- 
ing domain size suggests an obvious explanation for the 
appearance of the double transition upon annealing of 
as-grown samples (see Fig. [I]), this explanation appears 
to be inconsistent with magnetic x-ray and neutron scat- 
tering data which are interpreted as showing no change 
in domain size with annealing. 

In contrast, Model II is not as sensitive to the density of 
nanodomains. Specific heat signatures remain sharp even 
for a high density of 'nanodomains' and resemble those 
for high-quality crystals. Although these signatures re- 
main sharp, for coverages larger than ~ 75% only a single 
superconducting transition occurs. 



□ "as-grown" (Brison et al.) 
o annealed (Keizer et al.) 
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FIG. 1. Specific heat data of a low quality sample show- 
ing a single broad pealdp (squares) compared to that for 
a high-quality sample.E-3 The latter reveals signatures of 
two closely spaced phase transitions. In the normal state 
Cn/T w 430 mJ/(K 2 mol). 



While neither model provides an adequate account of 
the experimental data, the calculations presented here 
demonstrate the sensitivity of E 2ll superconductivity to 
SIAFM and suggest that a reasonably accurate descrip- 
tion of the spatial variation of the underlying AFM is 
crucial in drawing conclusions about the symmetry of 
the superconducting states. 



II. COUPLED SUPERCONDUCTIVITY AND 
ANTIFERROMAGNETISM 

A. The E-representation model 

We begin by reviewing the Ginzburg-Landau free en- 
ergy for an E 2 „ superconductor coupled to an SBF. We 
take the E 2 „ gap matrix A Q ^(k, r) to be the supercon- 
ducting order parameter; it has the form 



Aa/s(p/, r) = e(p/) ■ Tj(r)(c ■ iaa y ) a(i , 



(1) 



where e(p/) = (ei(p/), e 2 (p/)) are basis functions on 
the Fermi surface that transform among each other un- 
der the operations of the Dg/j symmetry group, rj(r) = 
(?7i(r), 7/ 2 (r)) are Cooper pair amplitudes that are func- 
tions of the Cooper pair center of mass r, c is a unit vector 
along the c-axis of the crystal, and <Xj are Pauli matrices 
in spin space. Explicit expressions for E 2u and other ba- 
sis functions may be found in, e.g., Refs. [7]|^ and |32|. The 
GL free energy density is 

/(r) = /bulk(r) + / grad (r) + /field (r). (2) 

Symmetry constrains the form of the bulk, gradicnt-and 
field contributions for an E^ u superconductor to beQ 



fb^ = a(T)\ri\ 2 +P 1 \r l \*+j32\vv\ 2 , 

/grad = Kl(A^)(A%)* + K2(D i rj i )(D j r)j)* 

+ K3(Dir)j)(Djr)i)* + K i {D z ri j )(D : ,ri j )* , 

/field 



s^lbl 2 



(3) 

(4) 
(5) 



where a(T) — cto(T — To), ao is a constant, To is the 
transition temperature, b = d x A is the magnetic field, 
and Dj = dj — i(2e/hc)Aj is the gauge-invariant gradi- 
ent. In the calculations that follow for a spatially varying 
SBF, we restrict ourselves to two dimensions and K4 plays 
no role. In weak-coupling BCS theory for an E 2u order 
parameter, the parameters K 2[ and K3 are small for Fermi 
surfaces with axial symmetry|3 we take k 2 = K3 = and 
write k = Ki for notational convenience. Weak-coupling 
BCS theory also predicts that /3 2 = /3i/2 independent of 
the shape of the Fermi surface and that the homogeneous 
equilibrium order parameter that minimizes the free en- 
ergy is doubly degenerate, breaks time-reversal symme- 
try, and is of the form r\ ~ (1, ±i). 

The orientation of the AFM order parameter N(r) 
may fluctuate dynamically. An estimate of the magnetic 
fluctuation time T mag , from the energy-resolution lim- 
ited magnetic Bragg peaks, obtained from elastic neutron 
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scattering experiments, gives r mag ~ 5OOps.0 Because 
this fluctuation time is slow compared to the characteris- 
tic time scale of the superconducting state t sc ~ h/A f+ 
50 ps, we take the SBF to arise from static AFM ordered 
and calculate equilibrium solutions of the GL functional 
in the presence of spatially varying AFM order. 

The transition temperature below which AFM order 
occurs is an order of magnitude larger than the super- 
conducting transition. In a mean-field description of the 
AFM order, the order parameter is well developed at 
the transition to superconductivity. Below the super- 
conducting transition, the neutron scattering data sug- 
gests that the modulus of the staggered magnetization 
decreases. Over the temperature range for which a GL 
theory of superconductivity is valid, the modulus of the 
staggered magnetization changes by less than 5%. Since 
this change is small, we take the staggered magnetiza- 
tion to be a fixed external field arid, neglect the effect of 
superconductivity on AFM orderO The leading contri- 
bution to the free energy from the coupling of AFM and 
superconductivity is second order in r\ and of the form 



sbf 



2ea MN-T? 



- cos 29 sin 29 
£aoT1 ' [ sin 29 cos 29 P 



(6) 



where N = (sin #(r),cos #(r)) is a normalized direction 
vector of the staggered magnetization, and e is a coupling 
constant proportional to the square of the staggered mag- 
netization. Symmetry also alio wg-a, coupling of the SBF 
to the gradient in the E2 U modcl,0'E3 so that / gla d is given 
by 



grad 



Ami 
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with 



ki(1±£ ± N 2 



(7) 



(8) 



This term, together with the c axis gradient terms, de- 
termines the kink and possible tetracritical point on the 
upper critical field curve. Because the magnitude of the 
gradient coupling to the SBF is small, being of the or- 
der e± ~ e ~ AT C /T C , it will not significantly affect our 
results, and so we neglect the (direct) coupling of N to 
rj through the gradient terms for the calculation of ther- 
modynamic properties in zero magnetic field. 

The symmetry breaking term in Eq. (||) is combined 
with Eq. (||) to give the free-energy density 

a_{T)\ m \ 2 +a + {T)\ V2 \ 2 + ea Q sin 29 (^ry* + c.c.) 



-K(A??i)(A%)* 



-1] 

8tt 



(9) 



with a±(T) = a(T) ± sa$ cos 29. For temperatures very 
near the normal-superconducting transition, the second 
order terms that include the coupling to the SBF dom- 
inate and 'real' phases of the form 77 ~ (cos# , s'm9)e l ' p , 
minimize the free energy. 



B. Single AFM Domain and Estimation of GL 
Parameters 



Without loss of generality, the salient features of cou- 
pling to uniform AFM (single infinite domain) can be 
seen by taking N(r) = (0, 1) and the coupling to the 
SBF to be positive, e > 0, so that / s bf favors 77 _L N. At 
the temperature T*° m = T + e, a phase transition occurs 
from the normal state to a spatially homogeneous super- 
conducting phase r\ cx (1,0). At a lower temperature 
T c llom = To — e/0, there is a second phase transition to a 
time-reversal symmetry breaking phase r\ cx (1, ±ir(T)), 
where = 0%l0i and r(T) is a function that .-grows 
rapidly and smoothly from to 1 as T is lowerednJ 

The specific heat jumps (per volume) at the two phase 
transitions T^° m and T c llom , measured relative to the nor- 
mal state, are calculated, from a derivative of the free 
energy and are given byE3 



AC ,hom =T hom Q 2/ 2/3i2 



AC 



horn 



T^ m ag/2 / 9 1 



(10) 
(11) 



with P12 = The ratio of the heat capacity jumps 



AC; 



horn 



(1 + 0) 



(12) 



Using these equations and thCpSnc^ific heat data, we esti- 
mate the ratio w 0.42-0.65Hli-J Noting that the weak- 
coupling value of lies comfortably in this range, we take 
= 1/2 in the calculations below. Further, the measured 
c axis H C 2 slopesQH of ~ 6.6 T/K imply within a coarse- 
grained analysis that n/ao « 50Knm 2 , i.e., a zero- 
temperature GL coherence length £ = y n/(aoT^° m ) w 
10 nm. This estimate of £ hi the basal plane is in very 
good agreement with other reports. c3 



III. SPATIALLY INHOMOGENEOUS 
ANTIFERROMAGNETISM 



We now consider the possibility that the orientation 
of N varies in the crystal lattice, N = N(r), and ex- 
plicitly investigate two cases: (1) abutting antiferromag- 
netic domains of uniform size, £ a fm, with orientations dis- 
tributed equally among the three possible q- vectors, and 
(2) randomly dispersed 'nanodomains' with characteris- 
tic dimensions of order the superconducting coherence 
length. For SI AFM, an analytic solution is generally no 
longer possible and it is necessary to solve the GL equa- 
tions numerically. The resulting superconducting state 
will be complicated, because of the competition between 
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the condensation energy gained by r\ orienting in direc- 
tions preferred by the SBF and the gradient energy cost 
to twist the orientation of the order parameter from do- 
main to domain. The response of an E2 U superconductor 
to these two models for SIAFM differs as described be- 
low. 



A. Preliminaries 



For numerical calculation, it is convenient to intro- 
duce a scaled order parameter f\ = rj/r/Q, where 770 is 
the modulus of the real phase solut ion of a homogeneous 
single-domain, r] = y aoet/2(3i2, that appears in the 
presence of uniform AFM order at a transition temper- 
ature T^° m = T + e. All temperatures are given in 
terms of a reduced temperature t = (T^ 111 —T)/e. Scal- 
ing Eq. (^|) to the magnitude of the free energy density 
of the homogeneous single-domain solution in the high- 
temperature phase, \fo(t)\ = (ao^) 2 / (4/3i2) , a dimen- 
sionless Ginzburg-Landau free energy density, /, is ob- 
tained which is of the form 

/ = -2[(1 -rsin 2 0)1*7! | 2 + (1 - rcos 2 0)fe| 2 

-rsin 29 Re - \\n\ A + ^fg^ x V*] 2 



+r(|A^i| 2 + |A*| a ), 



a |2 ' 



(13) 



where all lengths arc measured in terms of the SBF length 
£ e = ^/n/a a e, D l = ^ e D it and r = 2/f. 

For ease of calculation, we take a square computational 
mesh with a step size Air = Ay = 0.2 £ e . We tested our 
GL simulations on triangular lattices, which are more 
natural given the apparent hexagonal crystal symmetry 
of UPt3, and have found no qualitative differences for 
averaged quantities. Periodic boundary conditions were 
imposed. We take the GL parameters to be (3 = 1/2 and 
K = a e^. 

As the superconducting order parameter twists to 
accommodate the spatially inhomogencous symmetry 
breaking field, time-reversal symmetry breaking phases 
may appear in localized regions even for temperatures 
near the normal-superconducting phase boundary where 
'real' phases are expected to dominate. To detect the 
appearance of these phases, we calculate 



M orb (r) = ~ tj(r) x T7(r)* 



(14) 



which-.is (apart from a factor dependent on the gradient 
tcrma3) the spontaneous magnetization that arises from 
the internal orbital motion of a Cooper pair. Note that 
M or b is a real vector constrained by symmetry to point 
along ±c. For a homogeneous single domain, M or b van- 
ishes between the two transitions, as it must for a 'real' 
phase. At the lower transition, the second component of 



the order parameter begins to grow with a phase rela- 
tive to the first and M or b increases rapidly with decreas- 
ing temperature signaling broken time-reversal symme- 
try. As expected, the temperature dependence of M or b 
is consistent with a second order phase transition in mean 

field theory, |M orb (T)| ~ J T c h ° m - T, for T < T c llom . 

The calculation of free energies and spatially averaged 
quantities involves a summation over the lattice. We 
adopt the notation 



N N 



N~ 



(15) 



i=l j=l 



where N 2 is the number of computational mesh points 
and the position vector r = (x,y). In this notation, the 
spatially averaged orbital moment is (M or b(r)) P . In the 
case of Model I, we solved for several antiferromagncti- 
cally ordered domain sizes on a square lattice. For Model 
II we have found solutions of the GL equations for differ- 
ent sets of random configurations of the SBF and have 
not found any discernible differences for spatially aver- 
aged quantities. We attribute this 'self-averaging' of the 
quenched disorder of the randomized SBF to our rela- 
tively large system size (this would be exact for an infi- 
nite system) and to the large concentration of nanoscale 
defects (9% - 79%). 

quenched, afm ordered SBF 



FIG. 2. Structure of AFM domains assumed in Model I. 
For convenience a computational mesh with a square geome- 
try was chosen; explicit comparison with calculations on tri- 
angular meshes do not change our central conclusions. 

B. Model I - Checkerboard 

The spatially varying staggered moment configuration 
consists of abutting domains as shown in Fig. |^. The ori- 
entation pattern of N(r) was generated so that nearest- 
neighbor domains do not have an SBF of the same ori- 
entation and so that there is no net staggered magneti- 
zation. The influence of the SBF is greatest at high tem- 
perature and the nature of the superconducting phase 
depends on the size of the AFM domains relative to the 
superconducting coherence length. At sufficiently low 
temperature, the order parameter is rj(r) oc (l,i) to an 
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excellent approximation; for most purposes the coupling 
to the SBF leads to a negligible correction to the order 
parameter. 

The effect of domain size on the superconducting tran- 
sitions is apparent in Fig. |^. For large domains two 
reasonably sharp transitions are apparent. They smear 
rapidly as the domain size is decreased until for £ a f m ~ 
2 £ e only one transition appears for the size of our compu- 
tational mesh together with imposed periodic boundary 
conditions. A comparison of the numerical calculations 
of Fig. for our model in two dimensions with the spe- 
cific heat calculations for Garg's one-dimensional model 
presented in Fig. 2 of Ref. |l3| shows that these are ip-good 
agreement given the simplicity of the 'toy' model.E9 
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FIG. 3. Top: The specific heat for Model I in units of the 
upper specific heat jump of the spatially homogeneous sys- 
tem, AC+° m /T c h + m , for various AFM domain sizes measured 
in units of the SBF length £ e . The temperature t is measured 
relative to the upper superconducting transition and decreases 
in the positive x-direction, t = (T c h + m - T)/e. Bottom: The 
corresponding spatially averaged spontaneous magnetization 
(M or b(r)) r as a function of t and for the same domain sizes. 

The lower panel of Fig. ^ shows the temperature de- 
pendence of the spatially averaged spontaneous magne- 
tization which reflects the nature of the two phase tran- 
sitions. For an infinitely large single domain, (M or b(r)) r 
vanishes in the high temperature phase where the or- 
der parameter is 'real,' and rapidly increases below the 



lower phase transition, T^° m , as described above. As the 
domain size decreases, (M or b(r)) r becomes finite but re- 
mains small for temperatures below T c + and above T^° m . 
The time-reversal symmetry breaking phase appears with 
increasing strength as the order parameter tries to twist 
from domain to domain. For the largest domains, the 
specific heat shows two phase transitions (albeit rounded 
by spatial fluctuations), even though (M or b(r)) r is finite 
between the two transitions. For thp smallest domain 
sizes only one transition is apparento with an onset ap- 
proximately that of the transition temperature in the ab- 
sence of coupling to the SBF, Tq. 
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FIG. 4. Spatial cuts of the relative angle <j> between 
the order parameters r/i and 772 along the x axis for a fixed 
y-coordinate at temperatures t = 1 and t = 6. The AFM 
domain size is 10 £ e and (f> = Z (771, 772). The same parame- 
ters are used in Fig. |H| A relative angle of <j> — n/2 signals 
a superconducting phase that breaks time-reversal symmetry. 
The low-temperature phase (t = 6) breaks time-reversal sym- 
metry on average, while the high-temperature phase (t = 1) 
breaks time-reversal symmetry predominantly in the domain 
walls (indicated as vertical dashed lines). 

The twisting and flapping of the two superconducting 
order parameters 771 and 772 across the domain walls is ap- 
parent in the plots of the relative phase angle 4> between 



?/i and 772 in Figs. |i] and 
the relative phase in Fig. 



The flapping or unwinding of 
J follows the AFM on average, 
producing reasonable narrow (order ~ £ e ) superconduct- 
ing domain walls in registry with the AFM domain walls. 
Occasionally, the relative phase unwinds from to tt over 
an entire domain. Fig. |E| shows the spatially varying su- 
perconducting order parameter for the particular AFM 
configuration shown in Fig. || and a domain size of 10£ e . 
The spatial variation of the modulus of r/ shown at low 
and at high temperature in Fig. ^, panels (a) and (b), 
tracks the underlying AFM domain structure. As ex- 
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FIG. 5. Domain structure of the Model I superconducting state induced by the coupling to homogeneous AFM domains of 
size 10 £ e x 10 £ £ with three equivalent orientations, 0° , ±120° , of the SBF. Distances are in units of £ e . Panel (a) : Contour plot 
of \ri(r)\/(\rj(r)\) r at steps 0.999 (solid), 1.000 (dash), 1.001 (dot), 1.002 (dash-dot), and 1.003 (dash-dot-dot) at temperature 
t = 6. |fj| is minimum in the domain center. Panel (b): Contour plot of |r7(r)|/(|fj(r)|) r at steps 0.80 (solid), 0.85 (dash), 0.90 
(dot), 0.95 (dash-dot), 1.00 (dash-dot-dot), and 1.05 (light solid) at temperature t = 1. \t)\ is maximum in the domain center. 
Panels (c-d): The complex order parameter components 171 (light) and 772 (dark) are plotted as 2D vectors, where the relative 
size of the vector is proportional to its magnitude. In the low temperature phase, t — 6, in panel (c), 771 _L 772 on average, while 
in the high temperature phase in panel (d), 771H772 almost always. The orientation of the SBF is the same as in Fig. El 



pected, at high temperature, there is a suppression of 
superconductivity at the domain walls. At low temper- 
ature, there is a small reduction in the component of r] 
parallel to N in the ~ (1, i) phase preferred by the fourth 
order terms in the free energy and a negligibly small sup- 
pression of superconductivity occurs in the center of the 
domain. At high temperature, the components of the su- 



perconducting order parameter are 'real' in the interior 
of the domains. The orientation of the order parame- 
ter attempts to follow the AFM order on average. The 
orientation of r) does not follow that of N perfectly. A 
side-by-side comparison of Fig. |^(d) in comparison with 
Fig. ^ shows that even for these relatively large domains, 
the different orientation of the SBF in adjacent domains 
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forces a compromise in the orientation of rj as it twists 
from domain to domain. For example, r/ makes a 45° an- 
gle with respect to the SBF at the center of the domain 
(0 — 10,0 — 10) in Fig. ||(d). For temperatures below 
the second phase transition, the superconducting phase 
is essentially with a small perturbation from the 

AFM order. This is evident from Figs. |] and ||(c) which 
show a sizeable phase angle ~ n/2 between the rji and 772 
components. At high temperature, the order parameter 
is mildly suppressed at the domain walls where the rel- 
ative phase angle between the 771 and 772 components is 
sizeable, reaching ~ tt/2 at the corners, indicative of the 
appearance of time-reversal symmetry breaking phases 
in the domain walls. 




FIG. 6. The ratio of the contributions to the free en- 
ergy from the SBF coupling term and the gradient term, 
{/sbf(r)) r /{/grad(r))r, for Model I and for various domain 
sizes. Note that for the smallest domains the ratio is sig- 
nificantly larger at high temperatures, t < 3, reflecting a new 
(strongly disordered) state in which the twisting of the or- 
der parameter decouples from the AFM order and on average 
(H(r)) r ~(l,i). 

For various domain sizes, the nature of the spatially 
inhomogeneous superconducting state is reflected in the 
ratio (/ s bf(r))r/(/grad(r))r, where (/ s bf(r)) r measures, in 
an average way, the extent to which the superconduct- 
ing order parameter tracks the twisting of the SBF, and 
(/grad(r)) r measures the energy cost of twisting the or- 
der parameter. The ratio (/ s bf (r)) r /(/grad(r)) r plotted in 
Fig. ^| is largely temperature independent and insensitive 
to domain size near T c+1 i.e., £ — »■ 0, for large domains. 
For the smallest domain size, the ratio is significantly 
larger owing mostly to a smaller (/ gra d(r))r and the stiff- 
ness of the condensate. This suggests the appearance of 
a qualitatively diffcrcpt phase - a strongly disordered su- 
perconducting phaseO As can be seen from Fig. ||, near 
T c+ this phase differs from that for larger domains, be- 
cause time-reversal symmetry is broken globally and not 



just in the spatially restricted regions between supercon- 
ducting domains. 

Since GL theory does not predict the proper tempera- 
ture dependence of the specific heat, we apply a correc- 
tion to recover the correct temperature dependence near 
the double transitions so that GL results may be com- 
pared with experiment. We compute the specific heat 
sclf-consistently within weak-coupling BCS theory for a 
homogeneous single AFM domain for different sets of T c 
splittings in the clean limit and compare it with the GL 
result. A reasonable approximation to the BCS theory 
temperature dependence for a small temperature range 
below the normal-superconducting phase transition can 
be obtained by multiplying the homogeneous GL result 
by a factor proportional to et, as shown in Fig. 



2.0 r gl 

o BCS: T._/T ct =0.877, E=0.0428T„ 

a— a BCS: T./r' t =0.823, e=0.0595T„ 
o o BCS: T /T„=0.691, e=0.1 148T„ 




FIG. 7. The specific heat for an Ez n superconductor com- 
puted from weak-coupling BCS theory for various SBF cou- 
pling strengths e (symbols), and from GL theory for a single 
AFM domain. The BCS theory temperature dependence can 
be approximately recovered from the GL theory specific heat 
using a multiplicative factor that is linear in te. 

Our results arc shown together with specific heat ex- 
periments in Fig. H For ease of comparison, the experi- 
mental results have been scaled so that the peak in the 
specific heat at the lower phase transition agrees with the 
corresponding feature for a homogeneous SBF. Calcula- 
tions for domain sizes ~ 6£ e — 20£ E show two phase tran- 
sitions that are blurred by the spatial inhomogeneity of 
the magnetic order in a way that resembles experiments. 
From the comparison, we deduce that the coupling to the 
symmetry breaking field is small, e ~ 18 mK, compared 
to T c+ w 540 mK. 

The domain sizes consistent with specific heat experi- 
ments, 10£ e — 20£ e (~ 30£o — 60£o), are much larger than 
the uniformly sized 1£ £ — 2£ e (~ 3£o — 6£o) domains at- 
tributed to neutron scattering experiments. For the lat- 
ter small domains, only a single phase transition would 
appear in these calculations. This reflects the relative 
stiffness of the superconducting order parameter as com- 
pared to the condensation energy gained from orienting 
the superconducting order parameter in a direction fa- 



8 



vorablc to the local AFM order parameter. 
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FIG. 8. Calculated specific heat for Model I for a single 
domain and for various domain sizes in comparison with ex- 
periment (filled circles) normalized by Cn/Tl3 For ease of 
comparison, the experimental data scaled by a numerical fac- 
tor (open circles) is also presented. The T-dependence of the 
GL results has been corrected according to Fig. ^. 



C. Model II - 'Swiss Cheese' 

The superconducting state of Model I is very sensitive 
to the size of the AFM domains. To explore the response 
of an F,2u superconductor to different spatial AFM config- 
urations, we consider a model that represents what might 
be viewed as an opposite extreme from the abutting uni- 
formly sized domains of Model I. In Model II, nanoscale 
sized AFM domains permeate a large single ('infinite') 
AFM domain like the holes in swiss cheese. This pic- 
ture is motivated by the observation of intrinsic defects 
(most likely dislocation lines or stackirw-ifaults) in high 
quality and high purity UPt3 samplescHtil and by the ob- 
servation that neutron scattering data cannot distinguish 
between the commonly accepted picture of abutting uni- 
form domains with small moments and small domains 
with significantly larger moments. These intrinsic de- 
fects may act as nucleation centers for the random-field- 
likc symmetry breaking field, 'nanodomains,' and may 
provide a natural explanation for the linewidth broaden- 
ing of the AFM Bragg peak in reciprocal space as seen 
in neutron diffraction measurements. □ 

We used a standard pseudo-random number generator 
to uniformly distribute the nanodomains on our compu- 
tational mesh, and to assign a random orientation for 
the direction of N for each nanodomain relative to N 
of the 'background' AFM (0(r) = ±120°). On a mesh 
of 160 x 160 points (32£ e x 32£ e ) approximately 3400 
nanodomains (nanodefects) are needed to cover all three 
orientations of the SBF equally. An example configura- 
tion of the AFM order appears in Fig. ^. In Fig. |iy we 
show a typical distribution of nanodefects on a mesh of 
32£ e x 32£ e with a concentration of defects that covers 



approximately 44% of the mesh. Each nanoscale defect 
has a cross-shaped 5-point layout on the mesh. Inter- 
actions between adjacent defects or between clusters of 
defects were neglected. Fig. 11 shows the specific heat as 



a function of temperature obtained by heating from deep 
in the low temperature phase. Two phase transitions are 
signaled by heat capacity jumps that remain sharp, con- 
sistent with second order phase transitions, even in the 
presence of a high density of nanoscale defects. Increas- 
ing nanodomain density does lead to a reduction in the 
splitting of T c and only a single phase transition is ob- 
served for a nanodomain density larger than ~ 75%. The 
temperature evolution of the spatially averaged sponta- 
neous magnetization shown in Fig. [ll] indicates that the 
low temperature transition that separates two supercon- 
ducting states is second order. As expected, the rapid 
increase in the spontaneous magnetization is correlated 
with the appearance of a jump in the heat capacity, as 
shown in the top panel of Fig. [by. Note that time-reversal 
symmetry is broken for the single transition that occurs 
at high nanodomain densities. 

quenched, disordered random SBF 













A 




sis 











FIG. 9. The spatial orientation of the SBF for the ran- 
domly dispersed nanodomains (small arrows) in the presence 
of a uniform SBF in the background (big arrow) of Model II. 




FIG. 10. A typical spatial distribution of cross-shaped nan- 
odomains (nanodefects) covering approximately 44% of the 
32£ £ x 32£ e numerical mesh in the presence of a uniform AFM 
background. 
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FIG. 11. Top: The specific heat for Model II in units of 
AC7 hony T hom for the nanoscale defect mo del for various de- 
fect densities on a 32 £ e x 32 £ e lattice. Bottom: The orbital 
magnetization for the same model and set of parameters as 
in the top panel. 




FIG. 12. The ratio of the contributions to the free en- 
ergy from {/sbf(r))r and (/ gra d(r))r. A large ratio reflects 
a stiff condensate; it is energetically less favorable for the 
superconducting order parameter to twist near the nan- 
odomains. Note that the second transition can be described 
by t- ~ exp ([1 — x/x cr ] In 3), where x cr ~ 88% is the critical 
concentration of nanoscale defects, where there is only one 
superconducting transition. 



In contrast to Model I, the contribution to the free en- 
ergy from the symmetry breaking field dominates that 
from the gradient term as reflected in the large ratios 
(/sbf(r))r/(/grad(r)) r shown in Fig. ||. The supercon- 
ducting order parameter is on average aligned relative 
to the (background) AFM order except in a region ~ £ e 
around a nanodomain. 




T/T c+ 

FIG. 13. Calculated specific heat as a function of tempera- 
ture for Model II for various concentrations of nanodomains. 
shown in comparison with experiment (filled circles) normal- 
ized by Cn/tM For ease of comparison, the experimental 
data scaled by a numerical factor (open circles) is also pre- 
sented. The T-dependence of the GL result has been corrected 
according to Fig. [?]. 

In Fig. [II] we compare our results of the specific heat 
with measurements on high-quality crystals. Since the 
heat jumps are very insensitive to the density of nan- 
odomains, a wide range of e values and nanodomain con- 
centrations are consistent with experiments. In particu- 
lar, an SBF coupling e/To — 0.0633 requires a concentra- 
tion of roughly 36% in order to account for the observed 
sharp double transition in heat capacity measurements. 

It is natural to expect that disorder will drive the lower 
temperature phase transition from being second order to 
first order or possibly a glass transition. Our numerical 
heating and cooling cycles have shown that upon heating- 
up and crossing the low temperature phase transition, 
the entropy is always smooth and hence the transition 
is second order. However, when starting the cooling cy- 
cle above T c _ we find a glasslikc, frustrated and strongly 
disordered solution for the order parameter, which gives 
rise to a discontinuity in the entropy upon crossing T c _ . 
This discontinuity is consistent with a first order transi- 
tion. However, a comparison of the calculated free ener- 
gies shows that the latter is energetically less favorable 
than the solution with a smooth transition and signals 
that the glasslikc solution is metastable. Assuming that 
the metastable glasslike solution is experimentally ob- 
servable when rapidly cooling down, our calculations give 
a small latent heat £ = T C -AS = /iT c _AC_, where /j, 
is a numerical factor of order fj. ~ 0.01. In other words, 
the latent heat is a small fraction of the overall measured 



10 



specific heat, Q tts T C+ C N « 200mJ/mol, with i/Q < 1% 
or even less. In a carefully devised heat capacity mea- 
surement this small latent heat should be observable if 
indeed a glasslike phase transition occurs. 



IV. CONCLUSIONS 



We have explored the effect of spatially inhomogencous 
antifcrromagnetic order coupled to the superconducting 
order parameter of an ~&2u superconductor in two models 
representing limiting configurations of the AFM order 
parameter: (I) abutting AFM domains of uniform size 
equally distributed over the three possible orientations of 
the AFM order parameter, and (II) small domains with 
dimensions of the order of the superconducting coher- 
ence length (nanodomains), randomly dispersed through 
a single AFM domain. Our numerical solutions of the 
Ginzburg-Landau equations show that Model I is very 
sensitive to domain size. Phase transitions are rapidly 
broadened and smeared as the domain size is decreased. 
For domain sizes less than ~ 2£ e only one transition is ev- 
ident. The results of our calculations in two-dimensions 
are in qualitative agreement with those of the simple one- 
dimensional model of Garg. In contrast, Model II shows 
sharp phase transitions for all nanodomain densities up 
to the point where the double phase transition gives way 
to a single phase transition. Our calculations for Model 
II show that low-lying metastable states affect the ther- 
modynamic properties of the system as it is cooled, and 
suggest the possibility that the lower transition may be 
weakly first order. In contrast, our calculations for Model 
I show no evidence of a first order transition. Although 
both models can yield sharp phase transitions like those 
observed in the heat capacities of high-quality samples, 
neither of them can account for the change in the heat 
capacity on annealing if the magnetic moments and do- 
main sizes do not change as a result of annealing. The 
qualitative difference in the results for our two models 
does caution that simplistic models involving domains 
homogeneous in size do not rule out the possibility of 
E211 superconductivity in UPt3. The relative insensitiv- 
ity of the low temperature time-reversal symmetry break- 
ing phases to SIAFM in both models provides a natural 
explanation of how an E2 U superconductor can provide 
a good description of the gap structure of UPt3 at low 
temperature and therefore transport properties that are 
in good agreement with experiment. In the temperature 
region that includes the two phase transitions, Model I 
appears to be too sensitive to SIAFM while Model II is 
perhaps not sensitive enough. It is, thus, more likely that 
SIAFM is not arranged in abutting domains of uniform 
size, but rather, there is a distribution of AFM domain 
sizes peaked around a particular size. Further calcula- 
tions are required to properly consider this possibility, 
which we examine in forthcoming work. 
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